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Abstract 

The decoherence of quantum correlation is investigated in the Heisenberg spin system with the 
asymmetric anisotropic interactions. The quantum entanglement and discord are used to quantify 
the quantumness of the correlations. By the analytical and numerical methods, we find that 
quantum discord decays asymptotically in time under the effects of the independent local Markovian 
reservoirs. This is markedly different from the sudden change of the entanglement. Before the 
disappearance of the entanglement, the dynamic behaviour of quantum discord is very similar to 
that of the entanglement. It is also shown that the discord declines rapidly for the interacting spin 
system compared with the case of noninteracting qubits. At an arbitrary finite temperature, the 
nonzero thermal discord can be enhanced by the asymmetric anisotropic interactions which induce 
quantum fluctuations. 
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I. INTRODUCTION 



The quantum aspects of correlations in composite systems are a key issue in quantum 
information theory Quantum entanglement is extensively regarded as one useful mea- 
sure of nonlocal coherence. The inseparable (entangled) states are applicable in many parts 
of quantum information processing, like quantum teleportation, cryptography, dense coding 



and quantum computation 



2|-|4|. However, the recent theoretical and experimental devel- 



opments have discovered that some separable (nonentangled) mixed states are useful for 
improving performance in some tasks of quantum computers |5H7|. The entanglement is not 
the only type of quantum correlation. To generally quantify the quantumness of correla- 
tions contained in bipartite systems, Ollivier and Zurek have defined a measure known as 
n 

quantum discord |8J . It is largely accepted that the total correlations contained in quantum 
states can be described by quantum mutual information which is a sum of the quantum 
discord and classical correlation given by Henderson and Vedral 0, [loj]. Quantum discord 
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171 ]. By means of the thermodynamics 



12|. 



has received much attention in theoretical studies 

of correlations, quantum discord can exhibit signatures of quantum phase transitions 
As is well known, the composite system loses unavoidably the quantum coherence due to 
interacting with the environments 
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191 ] . The dynamics of the discord were also analyzed 



in the system of noninteracting qubits coupled to the environments [14|. The robustness of 
quantum discord is beneficial for the realization of scalable quantum computing in contrast 
to the fragility of the entanglement 20 1. 

Because of the fundamental and practical values, it is necessary to investigate the dynamic 
properties of quantum discord in some concrete systems of interacting spins. In semiconduc- 
tor quantum dots and single molecular magnets, the magnetic properties of these systems 
were usually studied by the Heisenberg spin models with the asymmetric anisotropic cou- 



plings such as Dzyaloshinskii-Moriya interactions |2lN25|. In this paper, the Markovian 
decoherence of the quantum discord is obtained by the master equation when the inde- 
pendent thermal reservoirs are considered. In section II, the time- dependent behaviour of 
quantum discord and entanglement are analytically and numerically analyzed. Compared 
with the case of noninteracting spin systems, the effects of spin interactions on the decoher- 
ence of quantum discord are taken into account. In section III, the thermal property of the 
quantum discord is also investigated. Finally, a simple discussion concludes the paper. 
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II. THE MARKOVIAN DYNAMICS OF QUANTUM CORRELATION 



To clearly describe the dynamics of quantum correlation, we use two types of measures 
including the entanglement of formation and quantum discord. For an arbitrary bipartite 
state pab-, the total correlations are expressed by quantum mutual information [261 ] 

Hpab) = Yl S ^ " S (Pab), (1) 

i=A,B 

where pi represents the reduced density matrix of subsystem i and S(p) = — Tr(plog 2 p) is 
the von Neumann entropy. Henderson and Vedral proposed one measure of bipartite classical 
correlation C(pab) based on a complete set of local projectors {n^} on one subsystem B. 
After the local measurements, the reduced state of subsystem A can be written by 

p k A = ^Tr B {(l A ® n k B ) PAB (i A ® n* )]. (2) 

Here p k is the measurement probability for the k—th local projector and 1a is the identity op- 
erator on subsytem A. When subsystem B is described in the two-dimensional Hilbert space 
{|0), |1)}, a complete set of the local measurements can be given by {Ii k B = V\k) (k\V^ , (k = 
0, 1)} where the arbitrary unitary operation V(9, 0) for 9 G [0, 7r] and G [0, 2tt] is 



cos 2 sm^e 



l p -i<t> 



v & <f>) = I ; > 2 j. (3) 

sin |e 9 — cos | 

Then the classical correlation in the bipartite quantum state pab can be given by 

C(pab) = S(p A ) - sup {S(pa\b)}, (4) 
{n|} 

where S(pa\b) = J2k=o iP k ^(PA.) * s the conditional entropy of subsystem A and 
sup{S(pa\b)} signifies the minimal value of the entropy with respect to a complete set 
of local measurements {n^}. Quantum discord is simply defined by 

Q(Pab) = I(pab) - C(pab). (5) 

In the following study, we draw on the master equation to analyze the time-dependent 
behaviour of quantum discord in the quantum system with the asymmetric anisotropic in- 
teraction. For a two-qubit system, the total Hamiltonian is written by H = H + Hj where 
the intrinsic Hamiltonian of two qubits H = ^(a A + <?%) anc ^ the interacting Hamiltonian 



Hi = fEa=i y z a A a B + D ■ ((Ta x &b)\- To lay stress on the effects of the interacting Hamil- 
tonian on quantum discord, we assume that H <C Hj. Here a a is the pauli operator. For 
the case of antiferromagnetic couplings J > 0, the asymmetric anisotropy D = De z is con- 
sidered. It is reasonable to assume that two qubits are coupled to its local thermal reservoirs 

n 

Ea,Eb [20]. Suppose that the initial state at t = is p(0) = pab(0) <8> (\0e a 0e b ){0e a Qe b |) 
where \0e a 0e b ) denotes the vacuum state of the two local reservoirs. The evolution of 
quantum state between qubit A and B is given by the master equation 

Pab (t) = -i[H, pab] + L(pab), (6) 

where the Lindbald operator 

Hp) = fa + l )^ 2a iP a i - p a i a i - a i a ip) 

i=A,B 

+ ^ili ( 2cr tP°7 - (X*i °t ~ °I °iP) ■ ( 7 ) 

Here = n is the mean number of the thermal reservoir and 7i = 7 signifies the rate of 
spontaneous emission for each qubit. In fact, the temperatures of the thermal reservoirs can 
be varied from zero values to small ones after the local operations. At very low temperatures, 
the initial state is approximately the ground state p(0). The mean number n = exp {ui E )k B T)-i 
could be a nonzero value where ue- 

We adopt the numerical approach to calculate the density matrix at any time. For 
numerical purposes, it is convenient to regard the components of density matrix p as being 
written as a single column vector p in Matlab. Then the pre-multiplication of the operator 
op can be written as spre(a)p where the function spre(a) can be given by the Kronecker 
product of the identity and matrix o. In the same way, the post-multiplication of the 
operator pa can also be written as spost(a)p. Therefore, the master equation in Eq. (6) 
can also be expressed by pABif) = £>Pab{Q) where the superoperator C is obtained by the 
pre-multiplication and post-multiplication of the operators. Because the superoperator C is 
time-independent and is small enough to be diagonalized numerically, the i — th component 
of p{t) can be calculated by pi(t) = J2j a ij ex P s jt- Here Sj is the eigenvalue of C and satisfies 
Cik = ^jUijSjUj^ . The coefficients = Uij^2 k U~^pk(0). The numerical calculation of 

and Sj is completed by the program of exact diagonalization of the superoperator C. 

From the point of view of the practice, the initial state Pab{®) between qubit A and B is 
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Pab (t) 



(9) 



usually chosen to be the ground state which is expanded in the product space of two qubits 

|^ = -L(|01) + e-|10)), (8) 

where e ia = ■ ^ i s clearly seen that the ground state is just one maximally entangled 

Bell state. The general expression of quantum state pABif) satisfies 

( u \ 

x z 
z* y 
\ v ) 

According to the result of [16], the analytical solution of quantum correlation can be given 
by 

Q(Pab) = S(p B ) - S(pab) + min {S , Si}. (10) 

The entropy is S = -(u\og 2 ^ + ylog 2 ^ + xlog 2 ^ + vlog 2 ^) and Si = 
~ J2 q =o i 1+( ~ 2 1 - >96 ' 1 log 2 1+ (~ 2 1 ) 9ei where 9i = a/(m + x — y — v) 2 + A\z\ 2 . Meanwhile, the en- 
tanglement of formation can be used to evaluate the quantum correlation, According to 
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29), the value of the entanglement is obtained by 

E(p AB ) = 2max{|2| - y/mj,0}. (11) 

Through the numerical approach, the decoherence of quantum correlation can be given by 
Fig. 1(a) and 1(b). If qubits A and B are in the ground state, the maximally entangled 
ground state has the maximal value of quantum correlation, E(pab) — Q(pab) = 1 at the 
time t = 0. In Fig. 1(a), it is clearly shown that the values of quantum discord Q(pab) are 
decreased and infinitely close to zero with time. The nonvanishing phenomenon of quantum 
discord is different from the sudden disappearance of the entanglement. In the vicinity of 
E(pab) — 0, the nonzero values of the discord can also quantify the nonclassical correlation 
for separable mixed states which make possible some tasks of quantum computers. However, 
it is noted that the decay of quantum discord is very similar to that of the entanglement 
before the time of the vanishing of entanglement. This aspect demonstrates that quantum 
discord can describe the nonlocal coherence like the entanglement. It is also valuable to 
investigate the case where qubits are initially in a separable state ^^(O) = |10). From Fig. 
1(b), we find that quantum discord and entanglement can be generated and then vanish in 
time. The behaviour arises from the coupling between qubit A and B. 



For a special case of J = 0, the evolution of quantum state pABit) can be de- 



scribed by a completely positive trace-preserving map [30(. For a general two-qubit 
state /?ab(0) = J2ki mn a mn,ki\k) A{ m \ ® I0s( n l; ^he evolved state in time can be writ- 
ten by PAsif) = Yukl,mn ^j=0 Ylf=0 a mn,M 

(K Aj \k} A (m\K j Aj ) ® (if Bi /|Z) B (n|iifi i ,) where the 



Kraus operators K l0 = J J^dO)^ + v/T^l^ll), K a = J^flO)^, K, 



^ i2 



T (vT^|0)i(0| + |l)i(l|) and K i3 = J^l)^. Here |0(l)) i is the ground (excited) 



2ra+lVV- L 1 I / * \ I / ^ ivA lo y 2n+l 

state of qubits i = A, B satisfying of |0(l))j = =f|0(1))j and p(t) = 1 — e 7t 2 + ' means the 
probability of the atom exchanging a quantum with the reservoir. If the initial state is \i/j) 9 , 
the expression of pABit) is similar to Eq. (9), 

U ~2n + r 2n + 1 ~ 2n + 1 ^ ~ 2n + P' 

_l r n(n + l)p 2 _ flp . _ (w_+l)p >. 
X "2 [ (2n + l) 2 +l 2n + l Jl 2n + 1 JJ ' 

z =^( 1 -p)»y = 1 - u - x - v - ( 12 ) 

Therefore, the angle 8i of the entropy S\ in Eq. (10) is analytically expressed by 9\ = 
(2n+i) 2 0- ~ p) 2 - F rom t ne above analysis, it is seen that quantum discord is independent 
on the asymmetric anisotropy D if the initial state is the ground state \ip g ). In Fig. 2, two 
cases of the decoherence of quantum discord are shown when the interactions between two 
qubits are J = and J = 1 respectively. It is found that the values of quantum discord 
rapidly decline to very small ones under the impact of the interaction J ^ 0. This manifests 
that the dramatic loss of nonclassical correlation happens when there is the strong interaction 
between qubits in the system. 



III. THERMAL DECOHERENCE OF QUANTUM DISCORD 

In general, states of qubits for quantum information processing are usually influenced 
by thermal temperatures. Under the effect of the temperatures, the thermal decoherence 
is unavoidable. Therefore, it is of value to study the properties of quantum discord for 
the thermal equilibrium state. In the spin model, the mixed state p(T) at any equilibrium 
temperature T is expressed by 

= \H IVWil expC-ej/fcF), (13) 

3 
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where the partition function Z = J2j exp(— ej/kT) and \if)j), ej are the j—th eigenstate of the 
Hamiltonian H and corresponding eigenvalue respectively. The exact diagonal expression of 
H is given by a set of eigenstates and corresponding eigenvalues, 

ei,2 =J(-i±M),|^ a > = -^(|01>±^|10»; 

e 3 =\ ~ oo, |V 3 ) = |00); e 4 = J - + u, |^> = |H>- (14) 

Here the parameter r) — 1 + %D. Because the thermal state has the same form in Eq. (9), 
the calculation of quantum discord is also solved by Eq. (10). The thermal properties of 
quantum discord are shown in Fig. 3. It is found that the asymmetric anisotropy D can 
increase the values of quantum discord at any temperature T. This is the reason that the 
asymmetric anisotropy along the z direction induces the quantum fluctuation in the XY 
plane which can enhance the quantum correlation. The decay of quantum discord occurs 
when the temperatures T are increased. For finite temperatures, the values of quantum 



discord are always nonzero. The result coincides with that of [31(. As is well known, the 
thermal entanglement is gradually decreased to zero with the increase of the temperatures 
281 ] . In the context, it is the fact that the robustness of quantum discord against the 
decoherence is helpful for the realization of quantum computing. 



IV. DISCUSSION 



Using quantum discord, we mainly investigate the decoherence of quantum correlation 
for two interacting qubits coupled to independent reservoirs. Compared with the case of 
noninteracting qubits, quantum correlation of the states in the interacting qubits system 
decays very rapidly. The values of quantum discord decrease asymptotically in time. This 
nonvanishing phenomenon is apparently different from the sudden change of the entangle- 
ment. It is found that the evolution behaviour of quantum discord is very similar to that of 
the entanglement before the sudden disappearance of the entanglement. This demonstrates 
that quantum discord for inseparable states can describe the nonlocal coherence like the 
entanglement. Moreover, after the sudden death of the entanglement, the nonzero values 
of quantum discord also manifest the existence of quantum correlation for separable mixed 
states. The thermal decoherence from the temperature is also considered. It is shown that 
the nonvanishing quantum discord at any finite temperature can be increased by the asym- 
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metric anisotropy. Quantum discord also declines asymptotically with the temperature. 
Therefore, quantum discord is a general quantifier for nonclassical correlation. 
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Figure caption 

Figure 1 

As two different quantifiers of quantum correlation, the quantum discord (solid line) 
and entanglement (dashed line) are numerically calculated and plotted as a function of the 
decoherence time. The mean number of the thermal reservoir is chosen to be n = 1 and the 
rate of spontaneous emission is 7 = 0.1. The parameters are chosen to be J = 1, D — 0.2 
and u = 0.1. (a). The qubits are initially in the ground state IV'ab(O)) = \ip g ); (b). The 
initial state of the qubits is a separable one IV'ab(O)) = |10). 

Figure 2 

The decoherence of quantum discord for the two-qubit system is illustrated under the 
independent reservoirs. The solid line denotes the case of two interacting qubits in the 
condition of J = 1, D — 0.2 and u = 0.1. The dashed line signifies the other case of two 
noninteracting qubits for J = 0. The mean number of the thermal reservoir is chosen to be 
n = 1 and the rate of spontaneous emission is 7 = 0.1. 

Figure 3 

The quantum discord at the thermal equilibrium state for J = 1 and 00 = 0.1 is shown 
with the changes of the temperature T and asymmetric anisotropic interaction D. 
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